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\^ • We consider the linearized nonsymmetric theory of gravitation (NGT) within the 

O , 

^D , background of an expanding universe and near a Schwarzschild metric. We show 

(N 

5-H ' that the theory always develops instabilities unless the linearized nonsymmetric la- 

O-i. 
■^ ' grangian reduces to a particular simple form. This theory contains a gauge invariant 



(N 



kinetic term, a mass term for the antisymmetric metric-field and a coupling with the 
Ricci curvature scalar. This form cannot be obtained within NGT. Next we discuss 
■^ ' NGT beyond linearized level and conjecture that the instabilities are not a relic of 

^D ■ the linearization, but are a general feature of the full theory. Finally we show that 

^^ ! one cannot add ad-hoc constraints to remove the instabilities as is possible with the 



instabilities found in NGT by Clayton. 
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bJ!). I. INTRODUCTION 

>'■ 

/\ ' So far, Einstein's general relativity (GR) has stood all direct experimental tests. The 

H : 

precession of the perihelium of Mercury, gravitational lensing and the redshift of light are 
measured in agreement with the theory at the percent level, just to mention a few. The crown 
on this experimental evidence, a direct measurement of gravitational waves, is expected 
within a few years \^. 

However there are also reasons to try to extend GR. For example the mysterious nature 
of dark energy and dark matter might become resolved within a modified theory of grav- 
ity. The standard example of a modified theory of gravity is Milgrom's Modified Newtonian 
Dynamics (MOND) [2], according to which Newton's law gets modified at very low acceler- 
ations, presenting an alternative to dark matter. MOND has recently been covariantized by 
Bekenstein J3|. 
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Another reason to try to extend GR is the notion of generahty. Within the framework of 
GR torsion is not included in a natural, geometric way. Within standard GR, any calculation 
of the connection (either by requiring metric compatibility, or by using the first order for- 
malism) leads to the (symmetric) Levi-Civita connection. Of course one is then free to add 
torsion, but there is no way that torsion follows naturally from the theory. An interesting 
generalization of GR would generate torsion in a purely geometric way, analogous to the way 
the Levi-Civita connection is generated in GR. 



ii 



The Nonsymmetric Gravitational Theory (NGT) |^ is an extension of GR that drops the 
standard axiom of GR that the metric is a symmetric tensor. Therefore we write the general, 
nonsymmetric metric as 

Qfj-u = Gfiu + Bfj_i^, (1) 

where G^^^ = g{fiu), B^^ = g^^i,] and (■) and [•] indicate normalized symmetrization and anti- 
symmetrization, respectively, as defined in Eq. ()A1|) . Indeed there is no physical principle 
that tells us that the metric should be symmetric and therefore such a generalization is very 
interesting to study. These kinds of theories where originally studied by Einstein p] p| in 
an attempt to unify GR with Maxwell's theory. However it turned out to be impossible to 
regain the Lorentz force law [7|] [8|] . The theory was revived by Moffat |9| , but now with the 
interpretation that the antisymmetric part of the metric really produces a new gravity-like 
force. This extra structure produces interesting results on the issues of dark energy and dark 
matter |li| [ll| |l2|. It will be clear that such a theory produces torsion in a very natural way. 
Furthermore note that a dynamical nonsymmetric metric field also naturally arises from non- 
linear cr-models of string theory |lj|. Unfortunately the nonsymmetric theory of gravitation 
suffers from all kinds of problems. The first of these is the notion of non-uniqueness. Because 
of the extra structure in NGT, there is no natural choice for the lagrangian. Without a 
guiding principle, the most general lagrangian that one can write down (with up to two 
time derivatives) has 9 extra terms apart from the usual Ricci scalar term (see Eq. ()A8|l ). 
Non-uniqueness also comes in another way: if we consider the expansion of the theory in 
terms of the i?-field, we could in general write 

Qfiu = Gfj,u + Bfj,i, + pB^j_aB'^y + aB G^u + 0{B ) 

where B^ = G^^^G^^B^yB^p and Ba^ = G^^B^p. p and a are undetermined parameters. 
The inverse is defined by g^'^Qnp = S'^. With these two extra parameters one sees that the 



linearized lagrangian is determined by 11 undetermined parameters coming from the full 
theory and the decomposition of the metric tensor (0) • 

Another problem with NGT, found by Damour, Deser and McCarthy ^], is the absence 
of gauge invariance 

B^u — *■ B^i, + d^Xu — <9^A^ (3) 

for the 5-field. This absence occurs, since in the linearized expansion one obtains, apart 
from the gauge invariant kinetic term 

— — if ; Haf3^ = daBp^i + dfsB^a + d-yBa/3 (4) 

also undifferentiated powers of B that couple to the background curvature. The absence of 
gauge invariance might lead to the propagation of ghosts, or unacceptable constraints on 
dynamical degrees of freedom. Fortunately these_problems can be relatively easily solved by 
the introduction of a mass term for the 5-field iM UM- 

n 

The final problem we discuss was found by Clayton [la|. One can show that, if one starts 
with a GR configuration and a small -B-field, the S-field will quickly grow and therefore 
the linearization of the 5-field does not make much sense. Since the field equations for the 
symmetric and antisymmetric part of the metric decouple at linearized level, Clayton argues 
that the instabilities might be a property of the linearized theory only. At higher orders, 
the 5-field backreacts on the symmetric metric and therefore the GR background should be 
considered as evolving. Such an evolving background leads to an increase of the degrees of 
freedom of the 5- field 17i| and these extra degrees of freedom then might stabilize the field. 

n 

Another (phenomenological) solution to these instabilities is given by Moffat [13|- The idea 
of this solution is to introduce an extra Lagrange multiplier in order that the unstable modes 
dynamically vanish. 

For the moment forgetting Clayton's result, in this paper we consider NGT linearized 
around a GR configuration. By explicitly constructing two different backgrounds (FLRW- 
universe in section IIIH and Schwarzschild metric in section IIV|) we show that the evolution 
of the i?-field is indeed unstable. By considering the most general form of the linearized 
lagrangian, we can explicitly point out which terms cause these instabilities. After having 
shown this we argue in section IV (^1 that instabilities are also present beyond linear level. 
Therefore our results can be extended to the full theory and should not be seen as a relic 
of the linearized theory. If these arguments are correct we must then conclude that even 
full, nonlinearized NGT suffers from instabilities. Finally we show in section IV Dl that the 



£ = V^ 



dynamical solution of Moffat to solve the instabilities found by Clayton does not solve the 
instabilities we find. 

In section|n]we briefly summarize the linearized theory, in sections IlIII and HVl we evaluate 
the evolution of the i?-field in two different GR backgrounds. In section |3 we analyze the 
results, discuss the theory beyond linear level and consider dynamically constrained NGT. 
In section IVII we draw our conclusions. In appendix ^ we state our notation and give the 
derivation of the linearized lagrangian. Appendix ^ summarizes the geometric quantities 
needed in the text. 

II. THE LINEARIZED LAGRANGIAN 

Since GR is very successful, it is natural to assume that any modification of the theory 
should be relatively small. Therefore we consider NGT in the limit of a small B, but an 
arbitrary G. The linearized Lagrangian of such a theory will in general have the following 
form (see Appendix IX|) : 

R + 2A- —H^ + (-m" + pR)B^ 
12 4 

n (5) 

-aR^,B^''B^''-^R^^,fsB^''B''^ +0{B^). 

Here the curvature terms R^au/3, R^iu and R all refer to the background, GR, curvature. The 
coefficient 7 is determined by the particular choice of the 'full' lagrangian, f3 is determined 
by the decomposition of the metric ((21), while a depends both on the full lagrangian and 
the decomposition. Naturally, different choices of these coefficients lead to different physical 
theories, so though in principle we can always make a decomposition of the metric in such a 
way that a and (3 vanish, we have no guiding principle that tells us they must be zero and 
therefore we will keep them arbitrary. The coefficient 7 cannot be set to zero within the first 
order (Palatini) formalism (see appendix El- Since the curvature couplings break the gauge 
invariance (jHI), the mass m^ has been added to prevent ghost modes to propagate |lJ|. This 
is not as artificial as it may sound. A mass is natural, since it is automatically generated in 
the presence of a nonzero cosmological constant (see Eq. ()A23|) ). We may assume that the 
mass term is generated by a cosmological constant and therefore we have today, 

m^ ~ A < 10'«^ GeV". (6) 

Note that this inequality is not necessarily true at all times, especially not when one considers 
an epoch of the early Universe, since the cosmological term may change during the evolution 



TABLE I: The scale factor and conformal time in different eras 



era 


a 


rj 


de Sitter inflation 


«- h]. 


v<-^ 


Radiation 


a = HjT] 


^ <V<Ve 


Matter 


a=||(r? + r?,)2 


ri>Ve 



of the Universe (for example during phase transitions). The field equations derived from the 
lagrangian (0) are 



- a{B''''R\ + B''^'R\) - 2-iB'''^R\''p + 0(5^) = 



(7) 



^^v „ ^^ ^J.v 



XG^, + 0{B^) = Q. 



We therefore see that to this order the field equations decouple and it makes sense to consider 
the symmetric background, to be just a GR background. The problem then reduces to the 
study of an antisymmetric tensor field coupled to GR. 



III. THE COSMOLOGICAL EVOLUTION OF THE B-FIELD 

In this section we study the behavior of the B-field in an expanding universe [l2| . Our 
background metric is given by the usual (conformal) Friedmann-Lemaitre-Robertson- Walker 
metric (FLRW): 

Gf,u = a{r]fr]^^, (8) 

where t]^,^ = diag(l, —1, —1, —1), f] is conformal time and a{ri) is the conformal scale factor. 
The conformal time is related to the standard cosmological time by, adi] = dt. The scale 
factor during the different cosmological eras is given in table HI where Hj ~ 10^^ GeV is the 
Hubble parameter during inflation and rje is the conformal time at matter-radiation equality. 
For the rest of this section it is important to note that Eq. © implies that 



Hi->m. 



(9) 



Of course, this bound can be amply satisfied even when the 5-field mass, m, is much greater 
than what is indicated in Eq. ©. 



A. The lagrangian and the field equations 



The linearized Lagrangian for the -B-field reads in conformal metric, 



(10) 



Here the geometric tensors, R, R^i,, R/^upa, derived from the metric (jHI) are given in appendix 
|B]and 7i is the conformal Hubble parameter given by 7i = a' /a, where a' = da/drj, denotes 
a derivative with respect to conformal time. This lagrangian leads to the following equations 
of motion: 

- 2aM'JZ\Bpx + n\Bxa) - ^a^i'JZ\''pBx. = . 
Next we act with rj^'^dx on (fTT|) divided by 2a^ to obtain 

a-2 ( - AaUn' + (12/5 + AaH") - (24/3 + 4a)H^ J E^o 

+ (- I2l3{n^ + W) + ^a'^') a-^r]''^dxB^p (12) 

- a(n\r]''Pd,B^x + 'R-\r]''%Bxo] - 2jr]''%(a-Y'"J^\./3Bx.) =0. 

The and the i components of this equation lead to the following two consistency conditions 
(2 = 1,2,3) 

Xri'WjBoj = (13) 

do{XB,o) + yv"'dkB,, = 0, (14) 

where we defined 

X = a^2 I (12/? + 2a)n'^ + (12/3 + 4a - 2-f)n' - Im^aA (15) 

y = a-^{12p + Aa- 2^)H^ + (12^ + 2a)n' - \m'^aA . (16) 

Next we consider the equations of motion (lllj). The Oi and the ij components give 

a'^r^P'^dp (d^Boi + doBi„ + d^B^o) + 2XB,o = (17) 



a-2 



r^^'^Sp d^Bij + d,B,„ + d,B,i ]-n[ doB,, + diBjo + d,B, 



+ 2yBi, = (1^ 



Equations (fT!?|l . (fT^ . (fT7j) and (fT^ determine the complete evolution of the B field. 



B. Solving for the 'electric' component 

For the rest of our discussion we will only consider the 'electric' component of the B 
field: Boi = Ei, where it is understood that E is a vector. We don't study the 'magnetic' 
component, since it turns out that it is completely regular and thus for our purposes not 
very interesting. If we assume A" 7^ 0, then we easily find, using (jT^. (jT^ and (fTTj) that 



7]''drdsE + dQ{y-^do{XE))\ - 2XE = 0. 



To study the behavior of this equation, we first rescale E: 

E 
to obtain 






dodo - ^S^'d.d, + Miff 
where the effective mass term is given by 



E = 0, 



^e// 



^ ^ 2y 43^2 



(19) 



(20) 



(21) 



(22) 



So we see from (PT|) that E behaves just as a massive vectorfield, as long as, y j X > 0. We 
see for example that in Fourier space we have 



^k' + Kf 



5o5o + ^k^ + M, 



Ep = 0. 



(23) 



where k is the momentum. Therefore, if y/X > 0, -E is a fluctuating, finite field and large 
values of k will be suppressed. On the other hand, if y/X < 0, for sufficiently large momenta, 
{y/X)k'^ + Meff < 0. For these momenta the field E grows exponentially without bounds. 
Such a situation is unstable and thus not physically viable. Therefore, if in such a case the 
transformation (|2m) is regular, the theory is unphysical. One could worry as well about the 
cases when, Mg^-j < 0. However on dimensional grounds, the effective mass squared scales 
in the worst case a,s,lhf. Such a scaling results in a standard power-law enhancement on 
super-Hubble scales jl9| and presents no problem. 

C. de Sitter inflation 

During de Sitter inflation (see table HI), Eq. (PT|) becomes (we use subscripts /, r or m for 
the flelds to indicate wether we are considering inflation, radiation or matter era) 



[dod^~5'^d,d,+MJ]Ei = Q, 



(24) 



where 

^7 = ^((47-12«-48/3) + |^). (25) 

So we see that during de Sitter inflation, the field is perfectly regular. Notice that any 
'negative mass' contribution indeed scales as 1/rf. 

During de Sitter infiation the rescaling ()20|) is nonsingular as long as 

2 
Tfl 

6a + 24/3-27-^^^0 (26) 

D. Radiation era 



During radiation era we obtain 



Er = 0. (27) 



Hfrn^T]* — 4(7 — a) 

Here Mr is the effective mass during radiation, which is given by a complicated expression. 
Fortunately its precise form is not important for us. It is sufficient to say that it indeed scales 
as 1/1]'^ and therefore it does not cause any problems. We see however, that we might have 
problems in this case with the sign of the coefficient in front of the spatial derivatives. For 
example if we look at the beginning of radiation era (77 = l/Hj) we see that if we want y/X 
to be positive, we need that im? / Hf is at least, 0{a — 7). In other words we approximately 
need: 

m > |a-7|if/ ~ |a-7| X 10^^ GeV, (28) 

which, unless |a — 7I is very small, contradicts Eq. (jH)). Therefore if we require y/X to be 
positive, we either need to fine-tune a or 7 such that a — 7 is sufficiently small to satisfy the 
bound (j^Hj) . or we could use the more natural requirement that = 7. The latter case will 
be examined in more detail in section UlI Fl 

For now we assume that ja — 7I is not unnaturally small, so we can effectively neglect the 
term, Hfirn?!]'^ in Eq ()27|1 at the beginning of radiation era and we get, y / X = — 1. Therefore 
we obtain an unstable equation of motion for E, implying that E grows exponentially during 
radiation era. Of course this growth is a problem only if the transformation (j2Uj) is regular. 
One readily checks that during radiation era, Eq. (j^UI) is non-singular if 

4(a - 7) + H'jm^T]'' ^ 0. (29) 

This is clearly the case, since we consider a — 7 to be dominant and nonzero. Therefore we 
conclude that, unless a — 7 is tiny, the equation of motion pTjl develops an instability at the 
beginning of radiation era. 



E. Matter era 



During matter era the equation of motion becomes 



a„a„ _ g?.^.. + .)'C^-^128(3. + 6/?-2-,)^,^g^ ^ ^^^, 



Em, = 0, 



(30) 



H]m'^{7]e + rifrf-* - 128(6/? + 7) 

where M^ is again a comphcated - but for our purposes not very interesting - expression. 
The factor in front of the spatial derivatives does not cause any problems, since the term 
i7|m^(?7e + ri)^T]~'^ dominates over the term containing the a, (3 and 7. The field rescaling 
(J2(J|) is during matter era non-singular if 



-768/3 - 1287 + Hjm 



2^2iVe + vY 



ri" 



7^0. 



(31) 



This relation always holds, since the last term dominates. 

F. Special cases 

The first special case we consider is a = 7. As was noted above, the case a = 7 is special, 
since now the instability at the beginning of radiation era appears to be solved. We notice 
that if a = 7 we have 



X = y = a-'[ (12/? + 27) (H' + W) 



a?w? 



(32) 



Now if we assume (j^ to hold (/3 7^ u^W^ ~ 2q;)), the only era where we have potentially 
problems is radiation. However during radiation era we find that 



7i„ + 7i„ = 



X = - 



m 



(33) 



and therefore the equation of motion (PTj) reduces to 



2„2 



dod^ - 5'^didj + ma 



E., = 



(a = 7). 



(34) 



So we find that the unrescaled field behaves perfectly regularly during radiation era. 

The next case we consider is the (local) situation where a = 7, A" — > and X' 7^ 0. This 
exactly corresponds to the cases where Eq. ()20|) becomes singular. Also notice that X -^ 
is only possible if /? > —7/6. In this case one finds that (fT4|l reduces to 



{12(3 + 2a){n" - n^)E = 



(35) 



10 

and thus E is constrained to become zero at one point in time. Such a constraint on a 
dynamical degree of freedom would normally be unacceptable. Here this is not the case 
however, since we also find from P^ that at that moment, 

Meff -^ oo. (36) 

Therefore E is dynamically pushed to zero (it becomes nondynamical) at the moment under 
consideration and the constraint ()35p presents no problem. 

The last special case we consider is the case where not only a = 7, but also (3 = —7/6. 
Now we find that 

x = y = -^ (37) 

and therefore during all eras the equation of motion reduces to the following simple form: 
Despite the elegant form, there is unfortunately nothing that tells us to constrain /3. 



E = o (« = 7 ; /5 = -|)- (38) 



IV. EVOLUTION IN A SCHWARZSCHILD BACKGROUND 

The other case we will consider is the behavior of the S-field in a Schwarzschild back- 
ground. Although the background is clearly spherically symmetric, there is no need for the 
S-field to also be spherically symmetric, so we will not explicitly make it so. The metric (in 
spherical coordinates) is given by 

G,. = diagjl - ^,-(1 - ^y\-r\-r' sin' e^ (39) 

where Mg is the Schwarzschild mass parameter. The equation of motion for the i?-field 
reduces in this spacetime to 

^dpiV^H'^'^n + ^^(^B^^ - 2^B''^R\^^ = 0, (40) 

with 



V^^ = r2sin^ (41) 

and the relevant components of the Riemann tensor are given in appendix El Notice that 
(HUj) is independent of a and (3 since the metric is Ricci-flat. We next act with d^ on (j^ 



and divide by y/—G to obtain the following consistency equation: 

2 

^d.B^^" - 2^d,{B''^R\''(,) 

2 4 2 (42) 

+ IIL^i- _ TLb'^''R\% + — cot OB''' - 27 cot eB'^PR\% = 
r r 2 
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The precise structure of these constraints is not very clear, so we consider the four differ- 
ent components separately. We obtain (the index a in Eqs. fl44Him]) indicates the angular 
directions: a = 2 = 6',a = 3 = = 3), 

TB^^ + 7^ cot ^S^o + SdiB^^ + 7^9a5'^° 

ncotOB^^ - SdoB^^ + ndaB^^ 

QB'^ - 7^9ofi'° - TZdiB^^ - Sd^B^^ 

QB^^ + S cot OB^'^ - ndoB'^^ - TZdiB^^ - Sd2B^^ 






(43) 





(44) 





(45) 


0, 


(46) 



where we have defined 



Q 


- r 1 r* 


n 


_ m2 2'yMs 
- 2 r-i 


s 


_ m^ 1 47M3 
- 2 1 r3 


T 


_ m^ i-yMs 



(47) 
(48) 
(49) 
(50) 



After a rescaling 



~ 7? ~ 7? 



S 5' 

we can combine the constraints in the following more elegant form: 



d^Bf"" + ^B^" + cot OB 



2u 



0. 



(51) 



(52) 



Next we consider the (unrescaled) equations of motion for the 5-field. Using ()4()j) we find 
that 



dpRp"" + -H^f"" + cot en'^'"' + 4:5[t.5'(. (s - n)B^^ + ASi'j^^ {s - n)B^^ + 27^5^^ = (53) 



[2 "3]^ 



We now focus on B'^^. After the rescaling (|5ip. we obtain from the 01-component of (j53|) : 



^ ( cot ed^ + d^d2 + d'^ds + 2s\ 5°^ 

+ cot eid'^B'^) + d\d2B^^ + dsB'^) 
+ cot 9{d^B^^) + d\d2B^^ + d3&^) = 

Next we use (|H^ to write the last two lines in terms of B^^ to obtain 



d'>do + d'di + ^ (d^d2 + cot 9d^ + d^d^^ + 27^ + 9^ (^) + ^d' 



B'" = 



(54) 



(55) 



12 
We rescale the field to get rid of tlie single derivative term. 

87MS + TTi^r'^ 
where we have defined 

(57) 



a/47 Ms — m?r^ 
Finally we divide (|3^ by (7°° and get 

(r - 2M,Y ^ ^ r - 2M, 7^ 



^2 r (S 



5o9o - ^ r^^^i-^i -^L' + M,^^^ 



S"^ = 0, (58) 



where the effective mass term is now given by 

2M. 



^e// = ( 1 - -r^ ) ( 27^ - (7"9i9iA ) . (59) 



and the operator 1? is defined to be 

L2 = 1 [d-,d2 + cot ^^2 + -^53^3) (60) 

When the solution of Eq. (jHHjl is written in a factorized form, B^^ = J2im^i"T.i^)'^im{(^,(p), 
the operator L^ generates a centrifugal barrier term, —l{l + l)/r^, where / = 0, 1, 2, .. is the 
multipole moment. Therefore we see that, depending on the sign of TZ/S, we have similar 
problems as with the cosmological solutions (PT|) . If TZ/S is positive (we only consider fields 
outside the Schwarzschild radius, r > 2Ms), we have a normal, well behaving field, since the 
field for high values of / is suppressed. On the other hand if TZ/S is negative, high values 
of / are no longer suppressed, and the field B^^ grows exponentially without a limit. Such a 
situation is clearly unstable. We evaluate TZ/S using ()48j) and ()49|1 and find that 

5 " mVS + 27M, ■ ^^^' 

We see that there are certainly finite values for r and Ms where this quantity becomes 

negative. In these regions the field is unstable, as was explained above. Note that 7Z/S is 
negative when, 

^ U^f fa' 7>0 

r <rcr= < ) \ 1 ■ (62) 

Since there is no reason to believe that 7 is particularly small and that m~^ should be less 
than the Schwarzschild radius, r^ = 2Ms, for all spherically symmetric mass distributions in 
the Universe, we conclude that the condition (J62|) will be satisfied at least somewhere in the 
Universe. 
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Since the B-Held should be stable for all values of Mg and r, for a sensible theory the sign 
of TZ/S should be independent of Mg and r. The only solution for this requirement is fixing 
7 = 0. In this case we have TZ/S = 1 and the equation of motion (|KH|l simplifies to 

O^Oo - '- =^^ ( ^^ =^^9i9i + L']+ M'f ) 5°i = (7 = 0), (63) 



(r- 


-2M,)/(r-2M,) \ \1 




r \ r ) ^VJ 




(r-2M.)/ 2r-4MA 




^^eff- ^ (^ 1 ^3 ) 



where 

r^ _ 9 A/T ^ / _ 9-r - 4 A/T \ 

(7 = 0). (64) 

Thus we see that (outside the Schwarzschild radius) the mass gets a positive enhancement. 
The field behaves perfectly regular and similar -but not equal- to a 'normal' massive vec- 
torfield in a Schwarzschild geometry. 

V. ANALYSIS 

A. Discussion of the evolution in a FRLW universe 

The calculations done in section lTTTl lead to the conclusion that, if we assume that the mass 
of the 5-field is much less then Hi and we want to avoid fine tuning, the couplings of the 
5-field to the background curvature, we must choose = 7. If we do not make this choice, 
then at the beginning of radiation era the B field becomes unstable and grows exponentially 
without a bound. We have also shown that the field equations are regular for this choice of 
the parameters. Looking back at the calculation we see that the reason these instabilities 
occur lies in the fact that there are two functions in front of the different components in the 
field equations (Af and 3^, see Eqs. (jl5flTH|) ). that could have a relative minus sign. This then 
leads to the possibility of the 'wrong' sign in front of the spatial derivatives ()58p . Indeed the 
choice a = 7 is the only possibility where X = y independent of the free parameters of the 
theory. 

The fact that we have a constraint on the couplings to the background, means that given 
a non-linearized lagrangian, one of the parameters of the decomposition is actually fixed. 
So, while we are not able to a priori tell how to decompose the metric into its symmetric 
and antisymmetric components, consistency of the theory during all epochs of the Universe 
gives us one constraint. Since there are two arbitrary parameters in (j2)) we are however not 
able to completely fix the decomposition. The special choice /3 = —7/6 for example would 
completely fix it, but except perhaps for simplicity, we are not aware of any physical reason 
why would this be the correct choice. 
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One could hope that looking at curved FLRW spacetimes might give a new constraint. 
These cases have also been checked and lead to essentially the same results. 

B. Discussion of the evolution in a Schwarzschild background 



In section IIVI we considered the evolution of a non-constrained -B-field in a spherically 
symmetric (Schwarzschild) background. We found that, unless we choose 7 = 0, the -B-field 
always has an unstable mode for particular values of Ms and r. Since there is no reason to 
exclude these values, we can only conclude that we should take 7 = for the theory to be 
consistent. The seeds for these inconsistencies lie, similarly as in the cosmological case, in 
the possibility of a relative minus sign between the coefficients of the different components 
in the field equations (in particular the coefficients TZ and S, see fHSHl^ ). Unfortunately, 
as we have shown in appendix 7 = is not possible within NGT in first order formalism. 
In fact one finds that an initial lagrangian that linearizes to a theory with 7 = also has 
the coefficient in front of H^ equal to zero. So this becomes a trivial, nondynamical theory. 

One might hope that perhaps a formulation of the theory in the second order formalism 
solves the problems and indeed one finds that now it is possible to find a full Lagrangian 
that gives 7 = 0. However in [l^ it is shown that in such a case, one always gets at some 
higher order n in B terms of the form RB"^. In the next section IV (yl we argue that these 
terms play exactly the same role as the RB^ terms we considered. 

Therefore the instabilities are present both in the first and second order formulations. We 
must thus conclude that there is no geometric full Lagrangian, as considered in appendix 1X1 
that correspond to a stable theory at linearized level. 

This result is no big surprise, since it was already noted by Clayton in [ly] that NGT has 
linearization instabilities. The nature of the instabilities we describe is however essentially 
different from the Cauchy instabilities Clayton finds. There are two proposed solutions 
for the instabilities Clayton discusses. First of all one might hope that corrections beyond 
linearized level might stabilize the theory. We address this issue in the next section where 
we are lead to the conjecture that our instabilities make NGT unstable beyond linearized 
level. Secondly, it was proposed in [lS| to modify the theory with lagrange multipliers that 
dynamically put the unstable modes to zero. In section IV Dl we show that these dynamical 
constraints do not resolve the instabilities we address. 
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C. Extending the results beyond linearized level: A conjecture 

Our results so far only apply to the linearized version of NGT. Since we showed that in this 
limit the field grows without a bound, one is lead to the conclusion that the linearization 
Ansatz of a small 5-field, is not anymore valid. We therefore need to consider NGT at 
higher orders. Already if we look at the quadratic order, we encounter many mathematical 
complexities, since (as is clear from (0)) Einstein's field equations get a modification of order 
B^. Therefore the standard OR solutions G^,y for the symmetric metric are not anymore a 
solution of the theory beyond linear level. Instead we assume we have (at quadratic order) 

gi,u) = G^, + f^,{B\H\G), (65) 

where fav is some unknown function of order B^ that vanishes if 5^ -^ 0. Therefore the 

n 

i?-field backreacts on the background metric. It has been argued by Clayton [ly] that such 
a backreaction might stabilize the theory, since the -B-field in such a case acquires additional 
degrees of freedom (similarly to 17[ ) ■ 

We will now look in a more detailed manner at the equations of motion for the S-field 
at higher order in the specific case where the GR background is the Schwarzschild metric. 
First of all notice that there will not be a quadratic contribution to the 5-field equations 
of motion. The reason is that there will be no cubic contributions to the lagrangian, just 
as there is no linear contribution. So therefore the first order we can consider is quartic in 
the lagrangian and therefore cubic in the field equations. Although the function f^y may be 
complicated, it is clear that a typical term that arises in these equations is {Rpaa/s refers to 
the GR solution), 

RpaapB^^'B^^^B'^r (66) 

An explicit calculation shows that this term is equal to 



2J^ ^ M.(2M-r) ^ M.(2M-r) ^._^, ^ 

_ (5i2)2^^^f (^13)2 Ms 2 Q _ ^B-'^f^MsT sin^ e] B^^r 

^ ' 2M. - r ^ ' 2M. - r v ) s . 



(67) 



Therefore we see that the different components of the B field appear with different coefficients 
in front of them. More importantly we see that, depending on the values of Ms and r, there 
could be a relative minus sign between various terms. Looking back at the calculations in 
section HV1 one sees that these relative minus signs are also present at linear level (although 
the coefficients differ, see Eq. (I5HJ) ). We also have seen that, depending on whether we have 
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such a relative minus sign, we may or may not get a relative minus sign between the Bq 
and the Of terms, and thus we may or may not get instabilities. The question now is if the 
relative minus signs in ()fj7|l implies the same for the cubic sector of the theory. 

First we assume we can solve the cubic equations of motion for only one of the components 
of the 5-field (which we call h) . Then, in analogy with section IIV[ we obtain an equation of 
motion of the following form (we assume that single derivative terms can be eliminated by 
a rescaling of h) 

(1 + aih'^)dlh + a2h + a^h^ + a^h{dohf = 0. (68) 

The precise value of the various a^'s depends on the specific form of the higher order la- 
grangian and the function f^i, in Eq. (fHKjl . 02 and 03 depend also on the momenta (spatial 
derivatives) of the field component. Eq. (jHTj) strongly suggests that the sign of the different 
a's is dependent on the values of Ms and r. Therefore we see that 02 and/or a^ can be 
negative. Since 1 0:2,3! can be as large as we want (because of the momentum dependence), 
we find by numerical analysis that ()68p always can develop instabilities (independent of the 
sign of «! or 04). Therefore, although we have not proven it, the similarity between (pTTj) 
and ()5H|1 indicates that, for certain choices of the parameters and coordinates, the necessary 
conditions for the instabilities to occur can always be met. Thus it appears that NGT is 
also unstable at cubic order. In fact since at higher orders in B there is no new structure to 
be expected, we are lead to the conjecture that NGT is unstable at all orders. 

The largest pitfall in the above argument is that it might not be possible to decouple the 
different modes of the i?-field. In such a case (jHH|) would change to an equation of the form 

(1 + aibl)d^bi + a2bi + a^b\ + aMd^hf = V[bi] , (69) 

where T> is some second order differential operator acting in general on all components 
hi 7^ 61 and with derivatives with respect to - possibly - all four coordinates. It is clear 
from the above discussion that the 'homogeneous' differential equation (with "D = 0) would 
give instabilities. The coupling between the differential equations induced by the operator 
T> in Eq. (jU^ may be such that the field components mutually stabilize each other. The 
question is whether this indeed can happen in a physically meaningful way. To investigate 
this consider the following toy model: 

fit) - kjm = ag'it) 

^ (70) 

9"{t)-klg{t)=l3f{t). 
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We expect this system to be representative of the dynamics of small S-field fluctuations 
around some general background metric which also includes nonlinear antisymmetric metric 
contributions. In Eq. (ffn|l a prime denotes a time derivative, kj and kg refer to the momenta 
of the field components / and g, a and f3 are functions of the background metric fields 
'G^jy' and ^B^jj\ possibly their first space or time derivatives (and thus they may be linear 
in kf and kg), and the parameters of the theory (M^ and r in the Schwarzschild case). In 
general the structure of the second spatial derivative terms in Eq. fl70p is not so simple 
(c/. Eq. (jiHl)), but the case considered here is representative for a more general situation. 
Indeed, the above discussion shows that the general second order spatial derivative terms 
reduce to the form ()7()|1 for certain ranges of the coordinates. Looking at ()54|1 shows that the 
coupling in Eq. (|7n|l is exactly the type of coupling we would expect between the different 
field components in NGT after we rescale away all single derivative terms. A more general 
coupling could also include pg"(t) + crgit), but this would not alter the following discussion. 
The system (j70|l gives rise to the following eigenvalue equation 

{co^ + kj){co^ + k^^)+u^aP = 0, (71) 

when the solution is decomposed as 

f{t) = X^OiC 

i=l 

4 

9{t) = "^biC 



4 

T (72) 



1=1 



Thus the system is unstable if at least for one root '^{uJi) < (3 indicates the imaginary 
part). We can solve (fTTj) to obtain 



\[-^-^-a,^^i^,^,.^-^^ 



.'.^-^-,^-^p^^^^,^,.,r-^^y ,T3) 

First note that if ^{u'^) ^ 0, then there is at least one root for which ^{oj) < is satisfied, 
and the system is unstable. So the necessary condition for the stability of the system is, 
53((x;^) = 0. From Eq. (f7S)) we then easily infer that this condition is obeyed provided. 



a 



or 

\2 



/3 > -{kf - kgf (74) 



a 



(3 < -{kf + kgY. (75) 



It is clear that it is impossible to satisfy the second equation for all values of the momenta. 
However the first equation can be satisfied if a(3 > 0. When '^{u'^) = is satisfied, the 
sufficient condition for stability is then, 3fJ(a;^) > (indicating the real part), which implies, 

kj + k^ + ap <0, {aP> 0) (76) 

Since we also require, aP > 0, both are impossible to satisfy simultaneously. We hence 
conclude that there always is an unstable root of Eq. ()7H|1 with '^{uj) < 0, which destabilizes 
the system (ffOjl . 

While the full analysis of cubic, coupled differential equations is much more difficult, 
based on the analysis of small fluctuations around a nonlinear background, we have argued 
that the same results apply: it is not possible to couple the unstable modes in such a way 
that the system stabilizes for all values of the momenta. The fact that the coupling constants 
in general depend on Mg and r only make matters worse. 

Finally, our analysis of the case where the different components of the 5-field mutually 
couple, indicates that our conjecture holds and NGT is indeed unstable beyond hnearized 
level. 

D. Dynamically constrained NGT 



In order to solve the instabilities discovered in NG' 
introduction of an extra term to the full Lagrangian 



? by Clayton [l6|, Moffat proposed the 
ig that looks like 



^DNGT = ^f^g^'^Ji^.cpu]- (77) 

Here J^ is some source vector and (pu plays the role of a Lagrange multiplier. This new term 
then leads to the dynamical constraint 



9 



'°] = 0. (7J 



Since these are exactly the unstable modes Clayton found, these instabilities are dynamically 
resolved. One might hope that this extra constraint also solves the instabilities we have 
found. 
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1. The cosmological case 

For the cosmological case we find from (fT3|) . (fT^ . (fTTj) and (fTHj) that (if ^ 7^ 0) the 
dynamics of the dynamically constrained theory are determined by 



^0^0 - 6'''drds - Hdo + 2a^y 



B,, = 0. 



After a rescaling 



one finds that 



\/a 



dodo - S'-'drds + Miff 



Bij = 0, 



which indeed has well behaving solutions. 



(79) 



(80) 



^1) 



2. The Schwarzschild case 



However for the Schwarzschild case we will now show that the problems remain, although 
in a somewhat different disguise. We focus on the 12-component of the field equations (pn|) 
and obtain 



d^do + d^ds + 27^ 



1 d23 , q2d31\ 



B'' + ds{d'B'-' + d'B 



0. 



(82) 



Using the unrescaled constraint equations (PH|) and (jlH|) we obtain 






12 I „ii 



B'^+g 



44)^-(^.| 



B^'+d'cotOB^' = 



(83) 
Next we rescale B^^ = f{r)B^'^ with some function /(r) that is determined by the require- 
ment that the single derivative terms in (jHHj) cancel: 



| + 5i|)/(r) + 2|9i/(r) 



0. 



^4) 



This equation can be solved, although the resulting expression is rather complicated. The 
final result for the equation of motion is 



(r-2M,)2 7e^^ r-2M, „ 
aoao ^ T^rdiOi L + mass term 



B^^ = 0, 



5) 



r2 S 

where the operator L^ is defined in (I60|l and the mass term is some complicated function 
of r. We clearly see that our troubles are not resolved, since once again the sign of TZ/S 
determines whether the solutions are stable. Exactly as in section llVl fcf. Eq (j61|) ). there are 
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certainly values of M^ and r where the TZ/S becomes negative. This results in an unstable 
solution for B^^. It can be shown that the function /(r) is also regular in this regime. 

As a final hope for stability one might try to impose even more Lagrange multipliers 
to also remove the unstable mode discussed above. However one finds that the mode B^^ 
exhibits the same behavior as B^^ and therefore should also be constrained to zero. Then 
one can show that the only remaining mode, B^^ is unstable exactly when B"^^ = B^^ = 0. 
Therefore we conclude that it is not possible to use dynamical constraints to remove the 
instabilities we found. 

VI. CONCLUSIONS AND OUTLOOK 

We consider the evolution of the Linearized NGT lagrangian in two different GR back- 
grounds. The first case we consider is a FLRW universe during three different cosmological 
eras. We show that, if the mass of the -B-field is much smaller then Hj, the field always 
undergoes an unbounded growth at the beginning of radiation era. The simplest solution for 
this instability problem is to fix the parameters a = 7 in the lagrangian ©. 

The second case we consider is a Schwarzschild background. Here the problems are much 
more severe. First of all the occurrence of instabilities does not depend on the smallness 
of the mass of the 5-field. Secondly the only way to remove the problems is by choosing 
7 = 0. Now this choice is not available in the first order formalism of NGT and therefore we 
conclude that the only stable linearized lagrangian must have the form 

C = y/^ln + 2A - —H^ + (-m^ + PR)b^] + 0{B'^) (86) 

and, as we show in Appendix ^ this form cannot be obtained by linearizing NGT. 

Since these results apply only for the linearized case, one might hope that considering 
higher order corrections stabilize the theory. One particular feature of higher order correc- 
tions is that Einstein's field equations get modified and therefore the standard GR solutions 
for the symmetric metric also get a modification of order B^. However, we have shown in 
Eq (jHTjl that the seeds for the instabilities are also present at higher order and therefore, if 
we can decouple the equations of motion for each mode, there is nothing that can prevent the 
instabilities from growing. If we cannot decouple the equations of motion, the situation is 
more subtle. Now in principle we could hope that two (or more) unstable fields could couple 
in such a way that they stabilize each other. We have showed however by considering a 
simple, but representative example (see Eq. (ffUjl ). that it is in general not possible to get rid 
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of the instabilities. While we do not know exactly the nonlinear field equations, this result is 
most probably applicable to any coupling between the field components one can construct, 
up to any order. Therefore, while our discussion does not explicitly prove that the theory is 
unstable beyond linearized level, we have gathered enough evidence to conjecture that NGT 
is unstable for all orders in B. 

One might also hope that modifying NGT could stabilize the theory. One modification 
of NGT that one can consider is dynamically constrained NGT [l8|. This version of the 
theory was introduced to solve the problems with linearization instabilities discovered by 

r~i 

Clayton [16] and therefore it might also resolve the instabilities we consider. We have 

shown that, while this modification does solve the instabilities in the FLRW background, 

the evolution in a Schwarzschild background remains unstable. We have also shown that we 

cannot add extra constraints in such a way that the unstable modes we have found disappear. 

Another, perhaps more promising, way to go is to try to find theories with a nonsymmetric 

metric that linearize to the form in Eq. ()86|) . One possibility might be to consider a theory 

on a hermitian, complex 4 dimensional manifold. In such a theory, the imaginary part of 

ihe metric must be antisymmetric and thus plays the role of the S-field. It is argued in 

20| that the extra symmetries on the complex manifold, could act as diffeomorphisms for 

the 5-field, thereby saving the gauge invariance Q- Such a theory has the properties we 

need for a physically viable theory, suggesting a possible direction to continue the study of 

nonsymmetric theories of gravity. 
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APPENDIX A: LINEARIZED NGT 

In this appendix we sketch the calculation of the linearized NGT lagrangian, starting 
from the most general, covariant, two derivative lagrangian. We work within the first order 
formalism, where the connection is considered to be a dynamical field and use the following 
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notation and conventions, 

g^i, : full, nonsymmetric metric 

ly" : full, nonsymmetric connection 

V : covariant derivative w.r.t. the full connection 

V : covariant derivative w.r.t. the zeroth order (Levi-Civita) connection 
Q : geometric tensors of the full theory 

R : geometric tensors of the zeroth order theory (GR) 

W, = iy[^„] 



1 



levrG 



N 



1 

(Al) 

and we use the following expansion of g^i, in terms of its symmetric (G) antisymmetric {B) 
part 

g^.. = G^, + B^, + pB^^B^, + aB^G^, + 0{B'') 

(A2) 

g^" = G"" + 5^^ + (1 - p)B^''BJ + aB^G^'" + 0{B^). 
This implies that 

V^g = V^{l + \{^--p + Aa)B'^ (A3) 

We use G to raise and lower indices. For other terms we use the subscript („) to indicate 
it is n^^ order in B (in general when there is no such subscript will mean a 'full' quantity). 
The Riemann tensor is defined in terms of the connection 

Q^«/3 = d^wi:^ - dpW!^^ - w:x^ + W'^pW^^^ (A4) 

and the Ricci tensor is also defined as usual 

Q,. = R\x. = d.w;:, - d^w;:, - w^,wx + w^^w;,. (as) 

Since the Riemann tensor ()A4|) does not have its usual symmetries (in fact it is only anti- 
symmetric in a and f3), we can make another, independent contraction (that would be zero 
in GR) 

P,. = R\,u = d,W^^ - d^W^^ (A6) 
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Finally we define our covariant derivatives by the 'H — ' relation: 



^i^Qap = di,9ap - W^^gpp - Wl^^g^p 



(A7) 



And we choose our Lagrangian in such a way that the other possible definitions of the 



connection do not give new information. Our starting lagrangian is 



l| 



C 



-99 



/.w 



(A^ 



+ dWpW, + 2A], 
where the parameters ai, 02, etc. are unconstrained J22] and A is the cosmological constant. 



1. Calculating the connection 



To calculate the connection we vary this lagrangian with respect to the connection and 
obtain 



-d„ 



^2 



^p"'^ 



i ( ^v9^ 



-99- ) + bid.iV^g^-^^) + 5'^d,iV^gn + 2a^S^^d, 
T^gl^-^^) + S^^d.iV^g^''-^)^ 



(v^^^'^^0 



+ v^{ - 9'''wp^ - 9"'w:;, + 9"'S'^w^, + g'^w^;, 






+ 



^p Act _|_ „pfi„(T\ \ I 9^„^M [per] 



ci g^'g-^'-gar, + 9"'9'''9va + '^C2g^'9""''9(c.^) 



+ I (9'' 9"^ 9a, + 9'"9^''9,a - 9^" 9'^ 9 a, - ^'^^^'^'^.a 



d{g^'^''^;-g^^'^^6'';)W, 



0. 



W.M 



(A9) 
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We contract this equation on p and t] to obtain 

A + 8ai - ^ + h^ dp{^/^gg^^''^) + v^j [h^ + 2ci - Cg - M)g^^''^W^ 

+ fea^t^'^V, - ^&3^'^^W^[H + '^C29^''9^'''^9iap)W^,.^ (AlO) 

and the contraction on a and t] gives 

SS.lv^^?^^'')) + (3 + 2ai + ^a2 - h^d^g^^"^^ + V^^Sg^'W^, + 2g^W, 

+ (361 + hj^g^^^'^WP^, + (26i + h,)g^P'^^W^ + ( - 62 + 2ci + cg + 3t/)(7('^^)W^, (All) 



f"^ ~2\ y^n^ y y ymj •'*' [^a] 

We only need to calculate the connection up to first order, since one can easily check that 
all second order contributions of the connection to the lagrangian ()A8|) can be written as a 
total derivative. From ()A10|) we find up to first order that 

W, = SVS^p (A12) 

where we have defined 

^ _ 2 + 8ai - f a2 + h ^^^^^ 



62 - 2(cl + c3) - 3d 

In principle one could think that we have special cases if either the nominator or the de- 
nominator of ()A13|) is zero. This means that (on shell) either W^ = 0, or V^B^p = 0. It 
turns out however, that in the linearized equations W^ and VB^p play exactly the same 
role. Therefore nothing special (in comparison to the case where S 7^ 0) happens in the 
cases where one of them is zero. A case that is special is when both the nominator and the 
denominator of ()A13|) are zero. However in this case we cannot proceed with the calculation 
of the connection. We therefore assume that this is not the case. We redefine the connection 
as follows 

W^^ = r^. - b'^pW^, (A14) 



(A15) 
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which means that Tf* , = 0, and use ()A10|) ) and ()A11|) to rewrite ()A9|) . Terms that are 
clearly 2nd order or higher have been dropped. 

+ h [g^'"'^rho[,,] + g^'^^^i) + 2(cl + c3)(7('^^)^(^-)^(„,)rf^,] 

+ hg^P'^WJ^ + Kl^S^g^P"^ - S^g^'^'^^W, + 0{B^) = 0, 
where we have defined 
^=[^ + 3«i) (^2 - 2(cl + c3) - 3rf) (^2 + 8ai - ^aa + 6i) 
K = (^2ai - iaa) (62 - 2(cl + c3) - 3d) (2 + 8ai - ^02 + 61) ^ - ^^2 + ^(ci + C3) + rf. 

Next we multiply ()A15|) with gp„ and substitute the resulting expression back in. For con- 
venience we also multiply our equation with gppgaa and obtain: 

- b2 [9{pp)^a + 9{aa)'['Jp^]) + 2(ci + C3)g(^x^)rf^c.] (Al6) 

The zeroth order solution of this equation clearly gives the well known Levi-Civita connection: 

''K^ = {;}• (A17) 

The first order equation can also be solved and gives 

/ \ (A18) 

where the coefficients in the last expression are given by 

(^ + i3)(6i + l) 



(\) = ^ 



A^ + AB- 2i32 
{A + ?>B){hi + l) 



A^ + AB- 2i32 
^ 2K + L (A19) 

^ = 2(1 - 62 + ci + C3) 
B = -2(ci + C3) 
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2. Calculating the linearized lagrangian 

Now that we now the connection, we can calculate the linearized version of the lagrangian 
()A8j) . It is clear that the zeroth order lagrangian is 

^^^C = ^/^{R + 2k). (A20) 

It is also not difficult to see that ()A12|1 implies that 

W£ = 0. (A21) 

for the quadratic order we find (all T's refer from now on to the first order connection) 

(2)£ = v/=g| (^ (^ - P + 2a)B'G^' + (1 - p)B^''bA R,, + \m^B^ 

+ (ci + C3)G^'G„^rf^,]rf,,j + {d- ^62 + \{ci + c,))w,w)j 
+ B^^ (^(1 + fei) VafJ,,] - v.rf;.,) + ( _ ^ - ^ai + ^a2 - ^fei) (9^W^, - d,W,) 

(A22) 
where the mass for the S-field is given by 

^'^' = A(^-p + 4a) (A23) 

Plugging in the expressions we found for the connection ()A18|) and removing a total derivative 
term results in 

(2)/: = v^^l {]^{\-p + 2a)B^G^'' + (1 - p)B^''bA R,, + ^m^E^ 

+ ($ + Ul)H^ + S(V^fi^,)(V"S,A)G^'^ + (^ - 2fi)(V'^S'"^)(V,5,.)| 

(A24) 
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where 



H. 



af3"/ 



$ 



</)2 + 2e0)(62 + ci + C3-i: 



E = d,T?e'' (1 - 62 + 2(cl + c3)) + S^ f rf - ^62 + ^(ci + C3) - ^L 

V o o o 

+ T.(2(e + h^e-^-- ^ai + ia2 - ^61 - |(ai + ^)) - ^(0 - 0(4(ci + C3) - 262 + 2) 
^ = (62-l)(e-0)' + 20(l + &i) 

(A25) 

van Nieuwenhuizen has proven a theorem 2l|, which states that in flat space the only 
consistent action for a massive antisymmetric tensor field is of the form 



'Bflat 



' 12 4 



(A26) 



If we want to make sure that ()A25|1 reduces to this form in flat space, the terms with the 
covariant derivatives have to be combined to a curvature tensor. This can be done with the 
following identity 



{y>^B-f'){VpB^,) + {VpB>^P){V''B^ 



fJ,(T J 



B^"'Ba''Rau - B^"'B''^R.^^B + total derivative 



(A27) 



With this identity we find that the requirement ()A26|) is satisfied if 

3 12 

S = ^ - 2fi 

and then the lagrangian reduces to its final form 

£ = V^ R + 2A- ^H^ + (jm^ + (3r)b^ 

- aR.^B^^'^BJ' - 7i?^,,^i?^''i?"'3] + 0(i?3), 

where 



(A28) 



(A29) 



a = p + H — 1 



(A30) 



7 = ^ 
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One can show that there is no choice for the parameters in the full lagrangian that has S = 
and simultaneously satisfies the constraints ()A28|) (what one finds is that if S = \E' — 21] = 0, 
then automatically $ + |fi = 0). Therefore we conclude that it is not possible to get rid of 
the 5-field coupling to the Riemann tensor. In other words, 7 = is not allowed within the 
theory. 

APPENDIX B: GEOMETRIC QUANTITIES 

We give an overview of the geometric quantities we use in the text. For the FRLW case 
we have 



{lo} = n, {?,} = m„ {i^=n6i 



Roo — —STi , 


Rq - 


4«' 


Rij = {W + 2n^)6i,, 


Ri = 


-\{n' + 2n')6i 


R = -^{n' + n'). 

(Ji 







(Bl) 



For the Schwarzschild case we have (the index a denotes 2,3): 

,2M, 



R'.\. = S}.6'- 



pi a _ ra rl ^s 
d2 3 _ X3 r2^^ 



(B2) 



The terms of the form R^^^i, (no summation over A) are symmetric in /i and z/ and therefore 
they do not contribute to our equations. 
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